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TURBULENT  WAKE  OF  A  SUBMERGED  FLAT  PLATE 


I.  INTRODUCTION 

In  order  to  characterize  turbulence  in  a  dynamic  situation  at  least  two 
local  scales  are  necessary,  a  turbulence  length  scale  and  a  turbulence 
velocity  scale.  This,  of  course,  does  not  uniquely  define  the  turbulent 
state  but  will  provide  sufficient  detail  to  characterize  the  flow  field 
averages:  average  velocity,  momentum  flux,  etc.  For  practical  calculations 
one  is  often  satisfied  with  a  description  of  the  average  motion.  This  is 
especially  true  when  we  consider  the  impractical! ty  of  determining  the 
instantaneous  fields  which  would  result  from  direct  simulation  of  the 
Navier-Stokes  equations  for  flow  fields  which  are  non-trivial.  For  problems 
which  involve  the  evolution  of  a  momentum  wake  in  the  vicinity  of  a  free 
surface,  the  computational  domain  is  extensive  and  we  are  therefore  required 
by  necessity  to  formulate  the  problem  in  such  a  way  as  to  make  any  numerical 
approach  practical.  That  is  to  say,  we  hopefully  seek  a  formulation  which 
is  "parabolic"  so  that  each  point  in  the  domain  need  be  calculated  only  once 
and  the  "solution"  can  be  generated  by  a  marching  procedure  through  the 
computational  domain.  We  would  like  to  eliminate  any  extensive  iteration. 

The  recent  work  of  Skop  (1),  which  is  based  upon  a  matched  asymptotic 

analysis  of  the  Navier-Stokes  equations  at  large  Reynolds  number,  provides 
Manuscript  approved  December  9,  1984. 
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us  with  a  conceptual  framework  in  which  to  imbed  our  formulation.  This  work 
has  established  the  parabolic  nature  of  the  momentum  wake  as  the  inner 
solution  to  the  problem  which  is  matched  to  the  well  known  outer  solution 
providing  the  Kelvin  wake.  In  consistency  with  this  we  seek  a  turbulence 
model  which  is  also  parabolic  and  which  can  be  conveniently  joined  to  the 
parabolic  averaged  momentum  equation  to  define  the  unknown  Reynolds 
stresses. 

While  there  is  no  guarantee  that  a  more  detailed  turbulence  model  will 
provide  a  better  description  than  a  simpler  one,  we  begin  with  the  assump¬ 
tion  that  at  least  a  two  equation  turbulence  model  is  required  for  this 

application.  This  is  based  on  the  fact  that  we  need  to  describe  at  least 
the  local  turbulence  velocity  and  local  turbulence  length  scale  as  indicated 

earlier.  The  choice  of  a  turbulence  model  Is  dictated  by  a  number  of 

factors  Including  ease  and  cost  of  Implementation  as  well  as  scientific 
issues  involving  the  physics  and  accuracy.  The  latter  Issues  are  difficult 
to  assess  because  turbulence  is  so  ill  understood.  One  can  hope  to  shed 
light  on  that  question  only  after  numerical  predictions  are  compared  with 
experiment  and  often  only  for  the  specific  applications  where  this  is  done. 

Turbulence  modeling  is  plagued  by  the  truncation  problem.  The  effects 
of  turbulence  upon  an  average  flow  field  is  at  the  most  fundamental  level 


due  to  the  interaction  of  the  Reynolds  stresses,  -  pu^u^  (second-order 
correlation)  in  the  averaged  momentum  equation.  The  evolution  equation  for 


the  Reynolds  stress  involves  third  order  correlations  such  as  u^u^u^  and 
other  unknown  correlations  like  pu,  ,  (pressure-strain  correlation).  Due  to 

1  >  j 

the  existence  of  the  nonlinear  inertial  term  in  the  Navier-Stokes  equations, 
evolution  equations  for  correlations  at  each  order  always  involve  higher 
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order  correlations •  This  approach  is  therefore  non-convergent  as  the  number 
of  unknown  correlations  always  exceeds  the  number  of  evolution  equations • 
This  divergent  procedure  must  be  truncated  at  some  level  otherwise  the 
system  becomes  hopelessly  complicated  without  any  guarantee  of  improved 
accuracy  or  an  ability  to  generate  solutions.  The  customary  procedure  is  to 
truncate  the  system  at  some  order  and  to  model  the  higher  order  correlations 
in  terms  of  lower  moments. 

Conventional  eddy  viscosity  formulations  fall  into  this  class.  Here, 
the  unknown  Reynolds  stresses  in  the  Navier-Stokes  equations  are  assumed  to 
be  proportional  to  the  mean  strain  rate  with  the  eddy  viscosity  the 
proportionality  factor.  The  turbulent  flow  is  thereby  modeled  as  an 
equivalent  laminar  flow.  The  eddy  viscosity  represents  an  enhanced  vis¬ 
cosity  to  simulate  the  stronger  turbulent  diffusion  processes.  The  eddy 
viscosity  so  defined  is  either  selected  to  be  some  representative  constant 
or  to  vary  throughout  the  flow  field  according  to  flow  field  geometry  and/or 
the  dynamics  of  the  turbulence  as  determined  from  some  modeled  turbulence 
evolution  equations.  One  approach  which  has  seen  wide  application  is  the 
two-equation  k-e  model  where  evolution  equations  for  k,  the  turbulence 
kinetic  energy,  and  e,  the  turbulence  dissipation  are  solved  along  with  the 
mean  momentum  equation.  The  eddy  viscosity  is  then  taken  as  some  function 
of  k  and  e.  The  most  common  form  is 


where  is  the  turbulent  eddy  viscosity  and  is  a  constant  typically 
taken  equal  to  0.09.  With  k  and  e  determined  throughout  the  flow 
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field,  v  is  generated  by  means  of  the  distribution  of  the  average 
turbulence  quantities  k  and  e  which  are  hopefully  in  balance  with  the 
dynamical  situation  in  question* 

The  earliest  applications  for  eddy  viscosity  were  ad  hoc  with  the 
existence  of  an  eddy  viscosity  assumed  as  an  ansatz  and  values 
of  v  prescribed.  More  formal  justifications  for  an  eddy  viscosity  and  the 
prescription  given  above  as  in  equation  (1.1)  will  be  established  later  in 
this  paper.  An  additional  complexity  exists  in  the  present  application, 
where  a  free  surface  forms  one  boundary  of  the  computational  domain.  It  is 
therefore  necessary  not  only  to  adequately  characterize  the  turbulence  in 
the  bulk  flow  but  also  to  characterize  the  turbulent  interaction  with  the 
free  surface.  We  must  recognize  the  possibiity  that  it  may  be  necessary  to 
model  the  turbulence  somewhat  differently  in  the  vicinity  of  the  free 
surface  boundary.  We  expect,  however,  that  at  remote  locations  relative  to 
the  free  surface  "standard"  turbulence  modeling  will  prevail.  The  free 
surface  interaction  will  require  free-surface  boundary  conditions  for  all 
computed  turbulence  quantities  and,  possibly,  proximity  effects.  That  is, 
there  could  exist  a  surface  layer  where  the  model  or  the  model  coefficients 
are  modified  from  those  in  the  bulk  flow. 

Very  little  work  is  available  on  the  problem  of  turbulent  flow  inter¬ 
action  with  free  surfaces.  That  work  which  is  available  seems  to  be  concen¬ 
trated  in  the  chemical  engineering  community  where  the  emphasis  is  on  mass 
transfer  through  interfaces,  i.e.  surface  renewal  (2),  (3).  Here  we  find  a 
number  of  visual  experiments  which  provide  some  insight  into  the  interaction 
of  turbulent  eddies  with  a  free  surface  and  accompanying  analyses  which 
support  the  conjecture  of  the  existence  of  an  interaction  layer. 
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Recently  Koraori  et  al  (5)  have  made  some  detailed  measurements  of  the 
normal  components  of  the  Reynolds  stress  in  the  immediate  proximity  of  a 
free  surface.  These  results  are  more  detailed  and  somewhat  contrast  earlier 
experiments  (4)  as  interpreted  by  Naot  and  Rodi  (6).  Current  evidence 
indicates  that  in  the  vicinity  of  a  free  surface  the  normal  fluctuations  are 
damped  while  the  lateral  fluctuations  are  enhanced.  The  enhancement  of 
lateral  fluctuations  is  consistent  with  conservation  of  turbulent  kinetic 
energy  at  the  interface.  Naot  and  Rodi  (6)  and  Hossain  (7)  have  introduced 
modifications  for  the  k-e  turbulence  model  by  means  of  an  algebraic  stress 
model  for  the  Reynolds  stress  with  free  surface  proximity  effects  similar  to 
those  developed  by  Shir  (8)  for  wall  proximity  and  a  free  surface  boundary 
condition  for  the  dissipation,  e  .  A  similar  model  was  used  by  Celik  et 
al  (9).  These  models  have  been  successfully  applied  to  open  channel  flow. 
Comparison  of  predictions  and  experiments  for  open  channel  flow  has  been 
good.  In  particular  the  expected  behavior  of  the  effective  eddy  viscosity 
near  the  free  surface  (14)  is  simulated  well.  Application  to  open  channels 
with  side  walls  correctly  simulates  the  secondary  motions  observed 
experimentally  and  predicts  the  observed  submergence  of  the  maximum  velocity 
point  below  the  free  surface  (10). 

While  open  channel  flows  differ  considerably  from  the  present 
application  of  turbulent  momentum  wakes,  these  are  to  the  writerfs  knowledge 
the  only  applications  of  turbulence  modeling  involving  free  surface 
effects.  Based  on  the  relative  success  of  these  former  applications,  in 
this  note  we  will  develop  this  particular  model  to  assess  its  applicability 
to  the  wake  decay  problem  in  the  vicinity  of  a  free  surface.  While  other 
turbulence  models  may  have  also  had  wide  application  and  comparison  with 
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experiment  for  different  applications,  the  inclusion  of  free  surface  effects 


seems  to  be  confined  to  Naot  and  Rodi  (6),  Hossain  (7)  and  Celik  et  al  (9), 
all  utilizing  similar  k-e  turbulence  models  with  an  algebraic  stress  model 
to  simulate  free  surface  proximity  effects. 

Although  differences  exist  throughout  the  literature,  turbulence  models 
are  usually  classified  in  terms  of  the  number  of  differential  equations 
which  are  utilized  to  describe  the  turbulence  quantitites.  An  extensive 
review  of  turbulence  modeling  has  been  provided  by  Rodi  (14).  From  what  has 
been  discussed  earlier  in  this  paper  we  consider  that  the  minimum  level  of 
description  for  an  arbitrary  flow  to  require  at  least  two  equations.  This 
is  because  of  the  need  to  describe  the  variation  of  a  turbulence  velocity 
scale  and  a  turbulence  length  scale.  It  should  be  emphasized  here  however 
that  the  particular  two  equations  need  not  characterize  the  turbulence 
velocity  and  length  scales  directly.  Any  two  turbulence  quantities  which 
can  define  the  local  velocity  and  length  scales  will  do  equally  as  well 
(14).  For  example,  the  k-e  set  provides  such  a  situation  where  the  velocity 


scale  u~  (-y*)  and  the  turbulence  length  scale  Z  ~  kJ/Z/e.  Recognizing 
this  together  with  the  fact  that  the  effective  diffusivity  for  momentum, 
etc.  scales  as  (11) 
v  ~  u  Z 


1/2 


3/2 


t 

we  recover  the  well  known  scaling  mentioned  earlier,  i.e. 
v  ~  k2  /e  . 

Historically,  zero-  and  one-equation  models  have  been  applied 
extensively  in  the  literature.  For  example,  constant  eddy  viscosity  models 
are  considered  to  be  zero-equation  models.  In  these  models,  a  constant 
value  for  the  eddy  viscosity  is  guessed  or  selected  by  a  consideration  of 
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the  overall  application  and  this  eddy  viscosity  is  used  throughout  the  flow 
field  to  relate  the  unknown  Reynolds  stresses  to  the  raean  velocity 
gradients.  The  "turbulent"  momentum  equation  is  solved  in  a  fashion  similar 
to  a  laminar  flow  field  with  the  constant  eddy  viscosity  used  throughout  to 
characterize  the  turbulence.  These  zero-equation  formulations  can  be 
expected  to  have  difficulties  in  situations  where  large  gradients  exist. 

For  those  applications  it  is  clear  that  average  turbulent  quantities  cannot 
accurately  describe  the  local  turbulence  values.  Improvement  is  to  be 
expected  by  the  use  of  a  one-equation  description  where  either  the  velocity 
scale  or  length  scale  variation  is  determined  through  a  differential 
equation.  For  applications  where  both  scales  vary,  both  variations  are 
necessary  to  provide  the  required  interaction  with  the  mean  flow  field  and 
we  can  anticipate  problems  with  one  equation  models  there  also. 

For  two-equation  models,  such  as  the  k-e  model,  there  are  two  classes 
of  approach.  In  the  earliest  applications,  the  local  values  of  k  and  e  were 
utilized  to  define  the  value  of  the  local  eddy  viscosity  by  means  of  an 
assumed  equation  of  the  form  (equation  (1.1)) 

v  =  C  k2/e . 
t  u 

Then  this  eddy  viscosity  was  used  to  relate  the  Reynolds  stresses  to  the 
local  mean  strain  rate.  This  then  provides  a  means  of  characterizing  a 
local  eddy  viscosity.  However,  it  has  been  established  that  many  turbulent 
flows  are  too  complex  to  be  described  by  a  simple  eddy  viscosity  concept,  in 
particular  flows  where  the  turbulence  itself  is  the  most  important  mechanism 
or  where  several  dynamical  effects  act  simultaneously  (12),  e.g.  shear  and 
buoyancy  and  possibly  for  our  application  where  free  surface  proximity 
exists.  To  circumvent  these  shortcomings,  more  recent  models  have  included 
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the  evolution  equation  for  the  Reynolds  stress.  This  equation,  of  course, 
contains  many  terms  and  usually  results  in  a  model  which  is  of  higher  order 
than  a  two-equation  model.  However,  the  inclusion  of  the  Reynolds  stress 
evolution  equation  provides  for  the  first  time  a  Reynolds  stress  relaxation 
time  scale  which  competes  with  the  hydrodynamic  time  scale  of  the  mean 
motion  (13). 

The  question  of  local  equilibrium  takes  on  a  more  meaningful  inter¬ 
pretation  and  further  provides  some  level  of  detail  into  such  questions  as 
isotropy.  In  terms  of  turbulence-free  surface  interaction  effects,  while 
various  experimenters  may  differ  in  their  detailed  measurements,  there 
appears  to  be  a  consensus  of  opinion  supporting  the  fact  that  the 
interaction  produces  an  anisotropic  distribution  of  turbulent  normal 
stresses  (4)  (5).  In  particular  those  normal  to  the  free  surface  are 
damped.  The  relative  enhancement  of  both  transverse  components  is  not  as 
universally  accepted . 

For  the  particular  application  which  we  are  stressing  here  it  would 
appear  that  some  consequences  of  the  Reynolds  stress  evolution  equation  may 
arise.  We  do  recognize  that  full  inclusion  of  this  equation  will  vastly 
increase  the  computational  requirements  and  will  additionally  require  models 
for  the  higher-order  turbulence  correlations  which  appear  in  this 
equation.  In  order  to  take  advantage  of  the  former  benefit  and  not  be 
penalized  by  the  latter  disadvantage  we  will  adopt  an  approximate  form  of 
the  equation  which  provides  an  algebraic  stress  model.  Therefore  no 
formidable  additional  computational  requirements  will  result  and  yet  we  will 
be  able  to  include  the  additional  physics  associated  with  a  free  surface 
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interaction,  i.e.  the  anisotropic  form  of  the  turbulence  normal  stresses  in 
the  vicinity  of  the  free  surface.  The  model  then  remains  a  two-equation 
model . 

It  seems,  therefore,  that  to  successfully  model  the  problem  of 
turbulent  flow  in  the  vicinity  of  a  free  surface  we  must  consider  several 
factors.  These  are  the  state  of  the  art  in  turbulence  modeling,  the  physics 
of  free  surface-turbulence  interactions,  the  necessity  of  maintaining 
control  over  the  required  computational  times  and  the  advantages  of 
employing  a  well-established  turbulence  model.  Then  we  are  led  in  a 
straightforward  manner  to  select  a  two-equation  k-e  turbulence  model  with  an 
algebraic  stress  model.  This  is  the  particular  type  of  formulation  that  has 
been  almost  exclusively  applied  to  model  free  surface  effects  in  open 
channel  flow.  Its  applicability  to  the  submerged  wake  problem  will  be 
addressed  below. 


II.  DEVELOPMENT  OF  THE  EQUATIONS 

In  this  section  we  will  define  the  concept  of  Reynolds  averaging  and 
then  develop  from  the  Navier  Stokes  equations  the  exact  averaged  equations 
which  it  provides.  These  equations  also  provide  a  starting  point  for  the 
generation  of  the  turbulence  model  equations.  To  simplify  our  discussion  we 
will  restrict  our  consideration  to  incompressible  flow.  This  assumption  is 
certainly  consistent  with  our  application.  For  this  case  the  Navier-S tokes 
equations  are  written  as 


3U1  „  3Ui  1  3  P  32“i 

3t  j  3XJ  p  3xt  3  XJ3  x  j 


(2.1) 
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The  continuity  equation 


completes  the  equation  set.  Here  represents  the  instantaneous  velocity 
component  in  the  x^  direction,  P  the  instantaneous  static  pressure, 
and  v  the  kinematic  viscosity.  For  a  turbulent  flow  application,  these 
equations  cannot  be  solved  exactly  for  a  practical  flow  situation.  In  order 
to  proceed,  a  statistical  approach  is  employed.  Toward  this  end,  the 
following  definitions  are  adopted.  All  dependent  variables  are  separated 
into  mean  and  fluctuating  components.  This  was  originally  suggested  by 
Osborne  Reynolds  many  years  ago  and  the  technique  to  be  described  below  is 
well  known  as  Reynolds  averaging.  We  therefore  take 


Ui  *  Ui+  ui  (2.3) 

P  *  P  +  p 

where  the  mean  quantities  and  P  are  defined  by  the  following  time 
averages 


Vh 


u 


i  dt 


(2.4) 


P 


t2  tl 


P  dt 


It  is  important  to  note  here  that  the  averaging  time  interval  t2-t^  has  to 
be  long  compared  with  the  time  scale  of  the  turbulent  fluctuations  and  small 
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compared  with  the  hydrodynamic  time  or  time  scale  of  the  mean  motion 
variations.  When  this  procedure  of  Reynolds  averaging  is  applied  to  the 
continuity  and  Navier-Stokes  equations  we  obtain  the  Reynolds-averaged  forms 


a  u  a  u 

— —  +  u  — - 
3e  J  to 


l  ap 

p 


a  x 


j 


3x, 


(v 


9  U, 


-  vv 


(2.5a) 


3x 


i 


0 


(2.5b) 


Equation  (2.5b)  also  implies 


9ui 

a  x , 


0.  Note  the  appearance  of  the  second- 


order  correlation  of  the  fluctuating  velocity  field  -  p  u^uj  the  momentum 
equation.  This  term  characterizes  the  turbulent  transport  of  momentum  and 
is  called  the  Reynolds  stress.  It  is  in  fact  the  appearance  of  this  unknown 
term  which  renders  the  equation  set  above  unsolvable  since  the  number  of 
unknowns  exceeds  the  number  of  equations.  To  generate  solutions  for  the 
above  equations,  this  term  must  be  modeled  or  expressed  in  terms  of  the 
other  variables  in  the  equation.  However,  we  can  generate  an  exact  equation 
for  the  second  order  correlation  by  means  of  the  following  procedure.  If  we 
subtract  the  time-averaged  momentum  equations  from  the  full  time-dependent 
Navier-Stokes  equation  an  equation  for  the  x^  component  of  the  fluctuating 
momentum  results.  If  this  equation  is  multiplied  by  the  fluctuating 
velocity  Uj  and  added  to  the  equation  resulting  from  the  same  procedure 
with  i  and  j  reversed,  and  the  equation  is  time  averaged,  we  obtain  the 


following  equation  describing  the  evolution  of  uiuj  (27) 
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3ul“j 


at 


+  u. 


3uiuj 


(ViuiJ 


-Ir!V 

P  3x^ 


3u  p 

+  -4-) 


3  x 


.1 


3  U, 


^iUk 


3  U, 


-  u  .  u.  5- 
.1  k  3 


*k 


3  u  3  u 

+  ^  [— 1  +  — i]  - 

P  3Xj  3  x£ 


1 2: 


(2.6) 


3  u  3  u  3  11  11 

2v  5-^  5— 1  +  v  » — ±-  J- 

3  xk  3  xk  3  xk3  xk 


We  note  the  appearance  in  this  equation  of  the  third-order  correlation 


u^UjU^  and  additional  terms  which  involve  the  correlation  of  pressure  and 
velocity  fluctuations.  For  high  Reynolds  number  applications,  the  last  term 


of  the  uiuj  equation  is  usually  neglected  and  henceforth  in  this  development 
it  will  be  omitted.  In  short,  our  ability  to  characterize  the  u^u^  equation 
has  not  really  improved  the  situation,  other  unknown  correlations  have  been 
introduced  which  have  to  be  characterized  before  a  problem  solution  can  be 
generated. 

It  is  the  task  of  turbulence  modeling  to  provide  a  tractable  algorithm 
from  this  situation.  The  basic  ingredients  of  turbulence  are  contained  in 
Equation  (2.6).  It  remains  to  provide  additional  information  which  includes 
the  relevant  physics  of  the  particular  application  in  sufficient  detail  to 
provide  a  legitimate  and  solveable  approximation.  Turbulence  modeling 
therefore  involves  the  art  of  approximation  while,  of  course,  balancing  the 
real  factors  of  computational  time  and  cost.  It  is  the  writer's  opinion, 
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shared  by  many  others,  that  a  "universal”  turbulence  model  does  not  yet 


exist.  Each  application  provides  for  its  own  optimum  model  reflecting  the 
economics  of  the  situation  as  well  as  the  basic  physical  complications. 


The  reader  should  note  that  upon  contraction,  the  u^u^  equation  becomes 
an  equation  for  the  turbulence  kinetic  energy  k,  where 


k  i  j  UjUj 


Therefore  the  turbulent  kinetic  energy  equation  is  obtained  as 


9  k 

97 


+ 


u.  9k 


k  3x, 


=_uiuk 


uiui, 


— jUk 


-  V 


(2.7) 


This  k  equation  or  a  variation  thereof  is  almost  universally  used  by 

turbulence  modelers  as  the  required  evolution  equation  for  the  turbulence 

velocity  scale  discussed  earlier.  To  complete  the  most  fundamental  two- 

equation  description  we  need  additionally  an  evolution  equation  for  the 

turbulent  length  scale.  While  some  variability  exists,  most  models  employ 

an  equation  for  the  turbulent  energy  dissipation,  e  (the  last  term  in 

equation  (2.7))  to  provide  this  local  information.  This  represents  a 

3/2 

legitimate  variable  since  the  scaling  suggests  e  ~  k  '  /l  .  Therefore  a 
combination  of  k-e  variables  enables  one  to  determine  the  local  turbulent 
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velocity  scale,  and  the  local  turbulent  length  scale  k^^/e  . 

Furthermore,  experience  has  shown  that  the  use  of  an  e  equation  simplifies 
some  of  the  modeling  problems  and  also  simplifies  the  considerations  near 
walls  (15). 

An  equation  must  therefore  be  developed  for  the  dissipation,  e .  It  is 
possible  to  derive  an  exact  equation  from  the  Navier-Stokes  equations  for 
the  fluctuating  vorticity  and  thereby  an  equation  for  e.  However,  this 
equation  contains  many  additional  complex  correlations  whose  behavior  would 
require  drastic  modeling  assumptions  to  be  of  any  use.  Therefore  the  model 
equation  for  e  is  usually  provided  somewhat  hueristically  rather  than  being 
derived  from  first  principles. 

The  modeled  form  of  the  e  equation  is  given  as  (14) 


3e 

3t 


+ 


—  3e 
U. 


i  3x, 


t  3e 

^  o  3  x 
e 


-) 

i 


+  C,  f  P  -  C, 
le  k  2e  k 


(2.8) 


Additionally  the  modeled  form  of  the  k  equation  is  presented  as 

3k  -  3k_  3 _  /_t  3k  % 

3t  i  3x^  3x^  0,  ^xi 


3  U 


3  U.  3  U 


Vtfe^  +  3^ 


(2.9) 
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2  _  3  U  3  U  3U.su  3  u  3  u 

with  v  »  C  —  ,  P  =  -  uu  x -  =  v  (* -  +  x — ^  x -  and  e  =  v  5 -  r — -  . 

t  \xe  i  j  3x.  tv3x,  3x,  3x,  3x, 


i 


j 


j 


1  1 


This  k-e  equation  set  is  seen  to  possess  five  constants  whose  specification 

becomes  part  of  the  model.  These  five  constants  are  C  ,  C,  ,  ,  a  and 

K  U  *  le  2e  e 

We  note  that  although  an  eddy  viscosity  type  formulation  has  been 
specified  to  model  the  diffusion  terms  in  the  k-e  model  equations,  the 
specification  of  the  Reynolds  stress  terms  in  the  momentum  equation  (2.5) 
has  not  yet  been  made  or  tied  to  an  eddy  viscosity.  Many  k-€  eddy  viscosity 
models  do,  in  fact,  employ  the  k  and  e  models  given  above  and  utilize  v  ^ 
given  above  as  an  eddy  viscosity  to  define  the  Reynolds  stresses.  This  is 
the  standard  k-e  model.  The  characterization  of  the  Reynolds  stresses  will 


in  general  arise  from  a  consideration  of  the  u^u^  equation  given  above. 
Additional  model  constants  may  be  introduced  at  this  stage  and,  most 
importantly,  additional  physics  may  be  introduced  as  well.  It  is  at  this 
level  where  we  will  introduce  the  experimentally-observed  phenomena  of  the 
anistropic  turbulence  Reynolds  stresses  in  the  vicinity  of  the  free  surface. 

There  have  been  various  attempts  recorded  in  the  literature  to  model 
the  Reynolds  stress  equation  and  include  it  in  a  turbulence  model.  However, 
it  should  be  emphasized  once  again  that  the  inclusion  of  this  equation 
within  a  turbulence  model  greatly  complicates  the  situation  and  would 
clearly  provide  the  need  for  more  extensive  computation  times  in  an  appli¬ 
cation  where  the  computational  domain  is  already  extensive.  Experience  has 
shown  that  higher  order  modeling  techniques  do  not  always  represent  an 
improvement  over  those  of  lower  order.  Our  consideration  of  the  Reynolds 
stress  equation  has  been  dictated  by  our  desire  to  incorporate  the 
turbulence-free  surface  interaction,  i.e.,  the  establishment  of  a  free 
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surface  proximity  effect  and/or  the  establishment  of  the  necessary  free 
surface  boundary  condition  required.  If  we  could  extract  this  information 
from  the  equation  without  including  the  full  differential  equation  as  part 
of  our  computational  load,  the  situation  would  be  ideal. 

Just  such  a  compromise  could  be  provided  by  an  algebraic  stress  model 
for  the  Reynolds  stress  equation.  This  approach  has  been  applied,  in  fact, 
to  open  channel  flow  in  a  number  of  investigations,  e.g.  Hossain  (7),  Naot 
and  Rodi  (6)  and  Celik  et  al  (9). 

In  Appendix  A  we  discuss  the  problems  of  modeling  the  Reynolds  stress 
equation.  Appendix  B  develops  the  required  algebraic  stress  equations  for 
the  Reynolds  stress  components  including  free  surface  proximity  effects. 

For  the  problem  under  consideration  here,  the  two-dimensional  submerged 
flat  plate  wake,  there  is  only  one  important  turbulent  shear  stress 
component,  uv«  This  component  is  approximated  by  equation  (B  21)  including 
free  surface  proximity  effects 


1 


3  u 


(2.10) 


where  f,  the  free  surface  proximity  term,  is  given  from  equation  (B.10) 


n  -  1  or  2 


(2.11) 
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and  l  is  a  local  dissipation  length  scale  of  the  turbulent  motion.  The 
geometry  is  illustrated  in  Figure  1. 

This  is  seen  to  represent  an  eddy  viscosity  formulation  for  the 
turbulent  shear  stress  -  uv  with  the  effective  eddy  viscosity  a  function  of 
distance  from  the  free  surface  y  through  the  proximity  term  f.  For  f  -►  0, 
at  remote  locations  from  the  free  surface,  we  recover  the  standard  eddy 
viscosity 


while  as  we  aproach  the  free  surface  y  +  0  the  effective  eddy  viscosity  is 
reduced.  This  is  consistant  with  open  channel  flow  measurements  (4).  For 
fixed  n  (linear  or  quadratic),  the  inclusion  of  free  surface  proximity 
effects  through  an  algebraic  stress  model  (Appendix  B)  is  seen  to  introduce 
two  additional  parameters  (CJ/C ^  and  Cf.  The  decay  of  the  eddy  viscosity 
near  the  free  surface  predicted  by  equation  (2.10)  seems  to  be  a  consequence 
of  modeling  the  anisotropic  normal  stress  behavior  through  a  Shir  (8)  type 
wall  proximity  function. 

The  complete  turbulence  model  description  is  seen  to  include  the  k 
and  e  model  equations  (2.8)  and  (2.9)  with  the  Reynolds  shear  stress 
component  given  by  equation  (2.10).  These  equations  are  to  be  solved  in 
conjunction  with  the  mean  flow  equations  (2.5a)  and  (2.5b).  As  is  typical 
with  high  Reynolds  number  turbulent  flow  in  the  absence  of  walls,  the 
viscous  stress  can  be  neglected  in  equation  (2.5a)  in  comparison  with  the 
turbulent  stress  (-  uv]  . 
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III.  BOUNDARY  CONDITIONS 


Boundary  conditions  are  necessary  for  all  the  dependent  variables  U,  k 
and  e.  These  conditions  are  required  at  an  "initial”  plane  or  upstream 
boundary  and  at  the  lateral  boundaries  of  the  domain  since  the  formulation 
is  parabolic  and  solutions  are  to  be  generated  by  a  marching  technique.  The 
lateral  boundaries  are  of  two  types  for  our  application  free  boundaries  and 
a  free  surface.  At  a  free  boundary,  free  stream  conditions  on  velocity  and 
turbulence  quantities  are  recovered.  Typically  this  implies  that  the  wake 
velocity  approaches  free  stream  and  k  and  e  vanish.  The  situation  at  the 
free  surface  is  not  so  clear. 

The  asymptotic  analysis  of  Skop  (1)  has  provided  the  proposition  that 
symmetry  conditions  are  to  be  satisfied  by  the  horizontal  velocity  at  the 
free  surface.  We  would  expect  further  that  the  turbulence  kinetic  energy 
would  also  satisfy  a  symmetry  condition  here.  The  experimental  evidence 
(4,5)  indicates  that  at  the  free  surface  there  occurs  a  redistribution  of 
the  normal  stresses  u2,  v2 ,  and  w2 ,  with  v2  damped  by  gravity  and  surface 
tension  and  with  u2  and  w2  enhanced  while  preserving  the  total  kinetic 
energy.  This  is  consistent  with  a  conservative  interaction  of  turbulent 
eddies  with  the  free  surface  or  a  symmetry  condition  for  k. 

The  boundary  condition  for  e  at  the  free  surface  is  more  complicated. 
Here  we  take  two  alternatives,  both  of  which  seem  to  have  some  experimental 
verification.  The  detailed  experiment  of  Koraori  et  al  (5)  and  their 
estimated  energy  budget  near  the  free  surface  suggest  that  £  has  zero  normal 
gradient  near  the  free  surface  for  open  channel  flows.  This  is  illustrated 
in  Fig.  2.  However  the  experiments  of  Ueda  et  al.  (4)  indicate  that  the 
eddy  viscosity  is  reduced  to  zero  near  the  free  surface.  This  seems  to 
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require  a  larger  dissipation  than  that  which  would  result  from  a  symmetry 
condition  for  e  and  suggests  the  imposition  of  an  enhanced  e  at  the  free 
surface.  The  experiments  of  Quarmby  and  Quirk  (24)  and  Hussain  and 
Reynolds  (28)  for  confined  flows  indicate  much  higher  measured  eddy 
viscosities  at  the  symmetry  plane. 

Hossain  (7),  Hossain  and  Rodi  (26),  and  later  Naot  and  Rodi  (16)  have 
suggested  that  while  the  dissipation  length  scale  is  reduced  as  we  approach 
a  free  surface,  similar  to  the  behavior  near  a  wall,  the  effective  dissipa¬ 
tion  length  does  not  go  to  zero  at  the  free  surface  but  goes  instead  to  some 
finite  limit,  £  .  Therefore  they  propose 


o 


where  k  is  the  turbulence  energy  at  the  surface  and  £  is  related  to  the 
s  ^  o 

open  channel  depth  D,  by 


£  *  0.07D. 
o 


(3.2) 


For  application  to  our  far  wake  problem,  we  anticipate  a  similar 
possibility  and  suggest  that  a  boundary  condition  for  e  be  given  as 


k 3/2 
s 

G  =  ~K - 

s  Ay 

m 


(3.3) 


where  ym  represents  the  depth  of  the  wake.  This  is  equivalent  to  defining 
the  effective  dissipation  length  at  the  free  surface  as 
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(3.4) 


l  ■  A  y 
o  m 

This  is  in  direct  analogy  with  equation  (3.2)  above  for  open  channels.  As  a 
free  surface  boundary  condition  for  e ,  two  alternatives  are  therefore 
suggested,  a  symmetry  condition  for  e  or  the  use  of  equation  (3.3).  Both 
conditions  will  be  utilized  in  the  calculations  to  be  presented  in  Section 
IV  and  their  plausibility  will  be  assessed. 


IV.  NUMERICAL  RESULTS  AND  COMPARISON  WITH  EXPERIMENT 

Above,  a  k-e  model  was  developed  to  describe  the  behavior  of  a 
turbulent  wake  in  the  vicinity  of  a  free  surface.  That  model  is  a  variant 
of  the  standard  k-e  model  with  the  inclusion  of  an  algebraic  stress  model  to 
describe  some  of  the  physics  of  turbulent  eddy-free  surface  interaction. 

For  two-dimensional  flow  that  model  is  summarized  by  standard  k  and  e 
transport  equations  given  by  (2.8)  and  (2.9)  and  an  eddy  viscosity  type 
formulation  for  the  Reynolds  stress  component,  uv ,  given  by  (2.10).  The 
eddy  viscosity  defined  by  (2.10)  possesses  surface  proximity  terms  f  as 
defined  by  equation  (2.11).  Without  the  surface  proximity  terms  or  in  the 
limit  as  f  +  0  we  obtain  the  standard  eddy  viscosity  relation 

v  =  C  k2  . 

Additionally  two  choices  are  provided  for  the  boundary  condition  for  e  to  be 
applied  at  the  free  surface,  e.g.  symmetry  or  £g  as  defined  from  equation 
(3*3).  Symmetry  boundary  conditions  are  taken  for  U  and  k  at  the  free 
surface.  The  model  can  therefore  be  viewed  as  possessing  three  levels  of 
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complication.  At  the  most  fundamental  level  we  have  a  standard  k-e  eddy 
diffusivity  model  with  symmetry  boundary  conditions  satisfied  by  all  the 
dependent  variables  at  the  free  surface.  At  the  second  level,  a  special 
boundary  condition  is  to  be  satisfied  by  e  at  the  free  surface  to  model  the 
reduction  of  turbulent  length  scale  there.  This  boundary  condition  has 
previously  been  successful  in  modeling  open  channel  flows  (2-6).  The  most 
complete  model  has,  in  addition  to  the  above,  a  surface  proximity  effect  to 
model  the  anisotropic  nature  of  the  normal  components  of  the  Reynolds  stress 
tensor  near  the  free  surface  (2-6)  as  has  been  observed  in  a  number  of 
experiments.  In  the  numerical  results  that  follow,  the  parabolic  forms  of 
the  momentum  and  k  -  e  equations  have  been  utilized.  This  is  consistent 
with  the  results  of  Skop  (1),  for  large  Reynolds  number  applications  as 
considered  here. 

In  this  section  we  will  attempt  to  optimize  the  above  model  in  terms  of 
the  proper  level  and  the  appropriate  free  parameters  by  means  of  a  compar¬ 
ison  of  the  numerical  predictions  that  it  provides  and  the  experimental 
turbulent  wake  data  of  Swean  and  Keramidas  (25).  They  measured  the  proper¬ 
ties  of  a  two-dimensional  turbulent  wake  induced  by  a  flat  plate  towed 
beneath  the  free  surface.  Measurements  were  made  at  several  locations 
downstream  of  the  plate  trailing  edge.  We  will  concentrate  on  the  first  two 
levels  of  the  model  as  described  above;  the  standard  k-e  model  with 
symmetric  boundary  conditions  and  the  special  boundary  condition  for  e  at 
the  free  surface  as  developed  for  open  channel  flows.  These  two  cases  will 
provide  insight  into  the  necessary  behavior  of  eddy  viscosity  in  the 
vicinity  of  the  free  surface.  Some  of  the  expected  consequences  of  the 
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surface  proximity  terms  can  be  determined  from  the  variation  of  the  e 
boundary  condition  as  both  tend  to  reduce  the  near-surface  eddy  viscosity. 

First  we  discuss  the  experimental  results  of  Swean  and  Keramidas. 
Second  we  present  the  results  of  model  parameter  variation  for  the  standard 
k-e  model  with  full  symmetry  conditions  at  the  free-surf ace .  Thirdly  we 
will  describe  the  model  results  associated  with  the  nonsymmetric  e  boundary 
condition.  From  these  results  on  the  e  boundary  condition  we  will  then  be 
able  to  establish  the  inapplicability  of  the  surface  proximity  effects  for 
the  two-dimensional  decaying  wake  problem. 

The  two  cases  above  represent  an  eddy  viscosity  formulation 
with  v  given  as 


For  applications  to  two-dimensional  flows,  there  is  only  one  important 
Reynolds  stress  component  and  it  therefore  is  reasonable  to  fix  and 

consider  only  the  variation  of  the  other  parameters.  In  addition  to  ,  the 
model  is  seen  to  contain  four  free  parameters  which  are  0£  ,  and  . 

a,  and  o  are  recognized  as  effective  turbulent  Prandtl  numbers  governing  the 

iC  £ 

diffusion  of  k  and  e  respectively.  C^and  are  model  constants  for  the 

dissipation  equation,  e.  For  the  purposes  of  this  paper,  we  will  fix  at 
its  universal  value  of  0.09  and  consider  only  the  variation  of  the  other 
four  parameters. 

It  should  be  noted  at  this  point  that  previous  applications  of  this 
standard  k-e  model  to  an  infinite  plane  wake  behind  a  flat  plate  by 
Rodi  (14)  have  provided  excellent  agreement  with  the  experiments  of  Chevray 
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and  Kovasznay  (30).  The  comparison  is  reproduced  in  Figure  3  to  provide 
some  confidence  in  the  standard  k-e  model  for  infinite  2-D  wake  flow 
computations.  The  values  for  the  parameters  of  Figure  3  are 


C 

V 


.09, 


°k  •  1-°  • 
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1.3  , 
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It 
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2e 


1.43, 

1.92. 


This  set  of  parameters  will  be  referred  to  as  the  "reference"  values.  They 
have  been  obtained  by  an  extensive  application  of  the  model  to  a  number  of 
types  of  free  turbulent  flows.  This  set  of  parameters  has  previously  also 
been  successfully  applied  to  many  two-dimensional  wall  boundary  layers,  duct 
flows,  free  shear  flows,  recirculating  flows  and  to  three-dimensional  wall 
boundary  layers,  confined  flows  and  jets  (14). 


A.  NRL  Experiments 

In  this  section  we  will  summarize  the  experimental  results  of  Swean  and 
Keramidas  (25)  as  these  provide  our  only  basis  for  comparison  with  model 
predictions.  The  plate  was  1  m  long  and  the  towing  speed  was  .59  ras~*  in 
all  cases.  Two  main  experimental  cases  were  investigated  for  the  two- 
dimensional  flat  plate  wake  flow.  The  first  case  which  is  called  "finite" 
corresponds  to  the  translating  plate  submerged  5  cm  below  the  free 
surface.  Data  were  collected  at  four  stations  downstream  of  the  plate 
trailing  edge.  In  the  second  case,  the  translating  plate  was  25  cm  below 
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the  undisturbed  free  surface  and  data  were  collected  at  the  same  four 
downstream  stations  of  the  trailing  edge.  This  case  is  called  the  infinite 
case  as  the  submergence  of  the  plate  was  sufficiently  large  so  that  the  wake 
flow  was  uninfluenced  by  the  presence  of  the  free  surface  for  the  extent  of 
experimental  domain  considered.  These  results  are  therefore  expected  to  be 
characteristic  of  an  infinite  domain  two-dimensional  flow  field. 

Figures  4-6  display  the  experimental  profiles  obtained  by  Swean  and 
Keramidas  (25)  for  the  infinite  case.  The  variation  of  mean  velocity 
relative  to  an  observer  fixed  on  the  plate,  turbulent  shear  stress  and 
turbulent  kinetic  energy  with  depth  measured  from  the  undisturbed  free 
surface  are  displayed  at  the  four  stations  downstream  of  the  trailing 
edge.  Figures  7-9  display  the  same  variations  as  measured  for  the  finite 
case.  An  obvious  difference  between  the  finite  and  infinite  cases  is 
immediately  apparent  from  an  examination  of  these  results.  The  infinite 
case  profiles  seem  to  display  a  permanent  anchor  at  of  the  depth  of 
submergence  of  the  plate.  That  is,  the  depths  of  the  maximum  velocity 
defect  point,  zero  shear  stress  point  and  the  central  wake  minimum  of 
turbulent  kinetic  energy  all  appear  to  be  fixed  at  the  original  plate 
submergence  depth  as  we  move  downstream.  Further,  these  profiles  are  nearly 
symmetric  about  the  plate  centerline  level.  This  is  not  the  case  for  the 
finite  experiments.  All  the  experimental  results  for  the  finite  case 
display  a  systematic  drift  of  these  above  points  toward  the  free  surface. 
That  is,  in  addition  to  the  expected  evolution  of  the  profile  form,  the 
entire  profile  appears  to  migrate  toward  the  free  surface.  The  data  at 
station  1 ,  x  »  5  cm.,  already  displays  a  "drift"  to  within  4.75  cm  of  the 
free  surface.  This  type  of  behavior  is  indicative  of  a  systematic  mean 
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vertical  motion  although  none  was  observed  experimentally.  Except  for  any 
local  elliptical  behavior  associated  with  the  free  surface  proximity  or  a 
drastically  asymmetric  vertical  diffusion,  this  migration  is  not  anticipated 
until  further  downstream. 

B.  Numerical  Results 

The  numerical  results  to  be  described  have  been  generated  by  means  of  a 
parabolic  Patankar-Spalding  (32)  code  as  provided  by  Rodi  (31).  This  code 
utilizes  a  forward  marching  finite  difference  procedure.  The  various 
turbulence  models  and  boundary  conditions  described  above  can  be  readily 
included. 

In  order  to  gain  a  direct  connection  with  the  experimental  data,  we 
utilize  as  upstream  boundary  conditions  profiles  obtained  from  the  actual 
experimental  data  at  the  first  experimental  station.  In  this  way,  we  hope 
to  eliminate  any  consideration  of  plate  free  surface  interaction  upstream 
and  immediately  downstream  of  the  trailing  edge.  We  will  then  be  presumably 
exercising  the  code  in  regions  where  it  is  most  likely  to  correspond  with 
parabolic  behavior.  That  is,  the  profiles  at  x  =  5  cm  have  been  obtained 
directly  from  the  measured  values  at  that  point.  Profiles  of  mean  velocity, 
turbulence  kinetic  energy  k  and  dissipation  e  are  thereby  specified  at  that 
station.  The  parabolic  Patankar-Spalding  code  is  then  utilized  to  propagate 
these  results  downstream  in  a  marching  fashion  to  produce  the  predicted 
profiles  at  the  other  three  experimental  stations  x  *  23,  62,  111  cm.  This 
facilitates  the  comparison  of  code  predictions  and  experimental 
observations. 
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1*  Standard  model-symmetric  boundary  conditions. 

In  this  section  we  will  describe  the  numberical  results  for  the 
standard  k-e  model  with  full  symmetry  boundary  conditions  as  applied  to  the 
two-dimensional,  submerged  flat  plate  problem.  We  will  provide  the 
predictions  for  the  two  experimental  cases  defined  above,  finite  and 
infinite.  In  both  situations  we  have  found  it  convenient  to  vary  four  of 
the  model  parameters  in  the  following  way:  (1)  for  fixed  values 
of  and  C2£  we  vary  the  effective  turbulent  Prandtl  numbers  and  for 

the  dif fusivities  of  k  and  e,  (2)  for  fixed  values  of  and  a  ^  we  then 

vary  parameters  of  the  £  equation,  and  .  Lastly,  (3),  we  consider 

the  effects  of  non-reference  values  for  all  four  parameters  in  a  few 
cases.  Our  main  objective  here  is  to  gain  an  appreciation  of  the 
sensitivity  of  the  results  and  structural  details  of  the  predicted  profiles 
to  the  various  model  parameters  so  that  an  optimum  parameter  set  may  be 
selected. 

i.  Infinite  Fluid 

The  results  obtained  for  the  infinite  fluid  case  will  now  be 
presented.  Figure  10  displays  for  the  reference  parameter  set  the  mean 
velocity  as  function  of  distance  below  the  undisturbed  free  surface  for  the 
four  x  -  stations  of  the  experiment.  The  reader  will  recall  that  the  x  =  5 
cm  station  represents  the  boundary  values  derived  from  experiment  which  are 
imposed  for  the  computation.  No  apparent  influence  of  the  free  surface  on 
the  results  is  obvious  here.  As  a  matter  of  fact,  there  is  no  significant 
fluid  deceleration  obvious  closer  than  approximately  18-19  cm  below  the 
free-surface  or  equivalently  beyond  6-7  cm  above  the  plate  centerline.  This 
is  further  exemplified  in  Figure  11  where  the  velocity  profiles  are 
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displayed  for  the  case  where  the  same  boundary  data  were  translated  to 
within  7  cm  of  the  free  surface.  The  subsequent  velocity  profiles  are 
essentially  identical.  This  translated  data  set  was  in  fact  the  one 
utilized  by  Swean  and  Keramidas  (25)  for  the  computations  in  their  paper. 
Figures  10  and  11  agree  reasonably  well  with  the  infinite  fluid 
experiments.  This  has  also  been  noted  by  Swean  and  Keramidas.  Therefore  we 
find  that  the  reference  parameters  do  provide  a  good  prediction  for  the 
velocity  field  data  in  the  infinite  case. 

Anticipating  some  of  the  results  to  be  discussed  later  for  the  finite 
fluid  case,  the  parameter  set  *  2.0  and  **  .75  with  taken  at 

their  reference  values  will  be  shown  below  to  provide  a  good  approximation 
to  the  finite  fluid  experimental  results.  Figure  12  shows  the  behavior  of 
the  velocity  field  for  this  parameter  set  as  applied  to  the  infinite  fluid 
boundary  data.  By  means  of  comparison  with  the  results  of  the  reference 
parameter  set  we  note  that  the  acceleration  of  the  minimum  velocity  is  only 
slightly  retarded  for  these  parameters  while  the  wake  structure  spread  in 
the  tails  of  the  velocity  profiles  appears  to  be  enhanced  for  this  case.  In 
the  infinite  fluid  case,  the  reference  parameters  appear  to  describe  the 
experimental  measurements  more  closely.  This  result  will  be  seen  to  be  at 
variance  with  our  application  to  the  finite  fluid  case.  This  factor  will  be 
discussed  in  more  detail  later. 

The  results  of  applying  the  same  parameter  set  as  in  Figure  12  but  with 
the  translated  boundary  data  set  of  Figure  11  is  displayed  in  Figure  13. 

For  these  parameters  we  note  there  is  free  surface  deceleration  (wake  drag) 
obvious  at  the  last  x-station,x  =  111  cm.  The  case  is  not  really  infinite 
for  this  translated  boundary  data  set.  Therefore  if  these  parameters  were 
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to  be  utilized  to  describe  an  infinite  case,  the  data  set  would  have  to  be 
translated  further  from  the  free  surface.  It  is  the  more  rapid  spread  of 
the  wake  tails  which  is  responsible  for  this  situation.  This  has  presumably 
occurred  because  of  our  utilization  of  an  enhanced  diffusivity  for  k  and  a 
retarded  diffusivity  for  e  relative  to  the  reference  case  in  our  modeling  of 
the  diffusion  terms  in  the  k  and  e  equation.  This  combination  is  equivalent 
to  an  enhanced  eddy  viscosity  away  from  the  plate  centerline. 

The  reference  parameter  set  seems  to  do  an  adequate  _1ob  in  describing 
the  wake  decay  for  the  infinite  case.  This  is  not  surprising  when  we 
consider  the  efforts  that  have  been  devoted  by  many  previous  investigations 
to  '’optimize"  this  reference  set.  By  definition,  the  infinite  case  should 
not  involve  any  free-surface  interaction  and  hence  the  reference  parameter 
set  which  has  been  obtained  from  the  examination  of  many  flows  which  do  not 
contain  free  surfaces  should  be  applicable. 

To  complete  the  description  for  the  infinite  fluid  case,  Figures  14-17 
display  the  evolution  of  the  turbulent  kinetic  energy,  k,  and  the  magnitude 
of  the  turbulent  shear  stress  for  two  of  the  parameter  sets;  the  reference 
set  and  the  set  which  is  found  to  be  satisfactory  for  the  finite  fluid  case 
to  be  considered  next. 

ii.  Finite  Fluid 

In  this  section,  we  will  consider  the  more  interesting  simulations  of 
the  code;  those  which  have  been  called  finite  fluid  since  the  plate 
submergence  is  sufficiently  small  relative  to  the  downstream  experimental 
domain  that  we  expect  a  significant  wake-free  surface  interaction.  As 
indicated  earlier  these  correspond  to  the  5  cm  submergence  case  measured  by 
Swean  and  Keramidas.  Once  again,  the  experimental  data  at  the  closest 
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downstream  station  is  utilized  to  define  the  upstream  boundary  values  of 
mean  velocity  turbulent  kinetic  energy,  k,  and  dissipation,  e .  The  boundary 
conditions  at  the  free  surface  are  taken  as  symmetry  conditions,  or 
vanishing  of  the  normal  gradients.  The  Patankar-Spalding  code 
with  k-e  turbulence  model  is  again  utilized  to  numerically  propagate  the 
solution  to  the  other  downstream  experimental  stations  where  a  comparison 
between  experimental  and  numerical  predictions  can  be  made. 

Figure  18  illustrates  the  results  for  the  velocity  field  obtained  for 
the  finite  boundary  data  propagated  with  a  turbulence  model  utilizing  the 
reference  parameters.  No  vertical  migration  of  the  maximum  velocity  defect 
point  is  predicted  within  the  downstream  domain  defined  by  the 
experiments.  This  will  be  seen  to  be  the  result  for  all  parameter  sets 
which  will  be  utilized  below.  It  can  be  seen  by  a  comparison  with  Figure  5 
that  this  reference  parameter  set  does  not  predict  the  free-surface  velocity 
very  well  and  seems  to  underestimate  the  value  of  maximum  velocity  defect  as 
compared  to  experiment  at  the  last  station,  x  a  111  cm. 

Figures  19-25  display  the  resulting  velocity  field  associated  with 
varying  o^and  while  keeping  C ^  and  at  their  reference  values.  It  is 

interesting  to  note  that  over  the  variation  of  and  considered,  the 
values  of  the  maximum  defect  velocity  are  essentially  constant.  It  appears 
that  these  parameters  preferentially  alter  the  tails  of  the  velocity 
distributions.  With  respect  to  variations  of  these  two  parameters  the  set 
associated  with  Figure  19,  =  2.0,  «  0.75  seems  to  describe  the 

experimental  trends  best. 

In  Figures  26-29,  and  are  returned  to  their  reference  values  and 
and  C2e  are  varied.  The  results  show  a  very  strong  sensitivity  to 
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variations  in  these  parameters.  The  effects  of  these  parameter  variations 
are  seen  to  modify  both  the  maximum  velocity  defect  and  free-surface 
velocity. 

Figures  30  and  31  illustrate  the  evolution  of  the  turbulent  kinetic 
energy,  k  and  shear  stress  profiles  for  the  reference  parameter  set. 

Figure  32  and  33  display  the  evolution  of  turbulent  kinetic  energy  and  shear 
stress  profiles  for  the  “optimum4*  finite  parameter  set.  This  set  seems  to 
model  the  experiment  more  closely. 

Table  I  summarizes  and  compares  for  all  parameter  sets  considered  the 
behavior  of  the  maximum  defect  point  and  free  surface  velocities  associated 
with  the  velocity  profile  evolution.  The  properties  of  the  experimental 
profiles  are  also  listed.  Parameter  set  2  seems  to  model  the  experimental 
profiles  best.  This  is  the  “optimum"  set  described  above. 

2.  Standard  model  -  e  boundary  condition. 

In  this  section  we  present  the  numerical  results  obtained  by  the 
application  of  the  standard  k-e  model  with  symmetry  conditions  imposed  at 
the  free  surface  for  all  dependent  variables  except  the  dissipation  e  .  For 
the  boundary  condition  on  e  at  the  free  surface,  equation  (3.3)  is  applied 
for  various  values  of  the  parameter  A.  This  boundary  condition  has  been 
applied  to  both  the  reference  parameter  set  and  the  “optimum"  set. 

Figures  34-43  illustrate  the  evolution  of  the  velocity  profiles, 
turbulent  kinetic  energy  distribution  and  turbulent  Reynolds  stress  as  a 
function  of  the  parameter  A  for  each  of  the  two  k-e  parameter  sets. 

The  results  display  for  the  smaller  values  of  A  an  enhanced  e  near  the 
free  surface  which  tends  to  reduce  k  at  the  free  surface  and  thereby  drive 
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TABLE  1 


COMPARISON  WITH  EXPERIMENT  OF  THE  EVOLUTION  OF  THE  PROPERTIES 
OF  THE  VELOCITY  PROFILES  FOR  VARIOUS  PARAMETER  SETS 


CASE 


(o 


k’ 


MODEL 

PARAMETERS 

a  ,  C,  ,  C„ 
e  ’  le  ’  2e 


MINIMUM  VELOCITY  FREE  SURFACE  VELOCITY 
@  (x2>  x3>  x^)  @  (x2>  x3,  x^) 


1 

(1.0,  1.287, 

1.43, 

1.92) 

(.67, 

.78, 

.832) 

(1.00, 

.996,  .972) 

2 

(.75,  2.00,  1 

.43,  1 

.92) 

(.67, 

.774, 

.83) 

(1.00, 

.986,  .96) 

3 

(1.00,  2.00, 

1.43, 

1.92) 

(.67, 

.776, 

.836) 

(1.00, 

.992,  .962) 

4 

(1.00,  1.00, 

1.43, 

1.92) 

(.67, 

.78, 

.834) 

(1.00, 

1.00,  .98) 

5 

(.75,  1.287, 

1.43, 

1.92) 

(.67, 

.77, 

.832) 

(1.00, 

.995,  .967) 

6 

(1.287,  1.287 

,  1.43 

,  1.92) 

(.67, 

.78, 

.832) 

(1.00, 

1.00,  .98) 

7 

(1.00,  1.50, 

1.43, 

1.92) 

(.67, 

.78, 

.832) 

(1.00, 

.995,  .972) 

8 

(1.00,  1.75, 

1.43, 

1.92) 

(.67, 

.778, 

.834) 

(1.00, 

.993,  .964) 

9 

(1.00,  1.287, 

1.43, 

1.75) 

(.66, 

.758, 

.81) 

(1.00, 

.998,  .987) 

10 

(1,00,  1.287, 

1.43, 

2.25) 

(.69, 

.81, 

.86) 

(1.00, 

.995,  .94) 

11 

(1.00,  1.287, 

1.00, 

1.92) 

(.70, 

.83, 

.88) 

(1.00, 

.976,  .922) 

12 

(1.00,  1.287, 

1.75, 

1.92) 

(.65, 

•74, 

.784) 

(1.00, 

1.00,  .993) 

Experiment  (25) 

(.66, 

.745, 

.82) 

(1.00, 

.99,  .95) 
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v  =*  k2/e  toward  zero.  Within  the  experimental  domain  this  effect  appears 
as  a  boundary  layer  concentrated  near  the  free  surface  which  affects  the 
surface  velocity  without  any  significant  interaction  with  the  maximum  defect 
points.  This  behavior  seems  to  reduce  the  predicted  deceleration  at  the 
free  surface.  To  provide  predictions  more  in  agreement  with  the  experiments 
large  values  of  A  ~  15  are  needed  which  actually  have  the  effect  of 
increasing  the  surface  eddy  viscosity.  This  behavior  is  at  variance  with 
that  provided  by  the  open  channel  experiments  where  best  agreement  was 
obtained  by  reducing  the  surface  eddy  viscosity  toward  zero.  The  decaying 
wake  seems  to  require  an  enhanced  eddy  viscosity. 

It  is  because  of  these  trends  with  the  e  boundary  condition  that  we  did 

not  pursue  the  results  for  the  free-surface  proximity  terms  tt  of  the 

1.1  *  s 

algebraic  stress  model.  The  effects  of  including  these  terms  in  the  eddy 
viscosity  should  be  similar  to  those  obtained  for  the  smaller  values  of  A 
since  both  reduce  the  eddy  viscosity  near  the  free  surface. 


V.  SUMMARY  AND  CONCLUSIONS 

A  general  two-equation  turbulence  model  has  been  developed  for 
application  to  wake  studies  in  the  vicinity  of  free  surface.  The  model 
includes  standard  k  and  e  equations  and  ability  to  include  surface  proximity 
effects  through  the  use  of  an  algebraic  stress  model  very  similar  to  that 
previously  applied  to  model  open  channel  flows.  The  entire  problem  is  cast 
in  a  parabolic  form  which  lends  itself  to  a  marching  procedure  in  which 
profiles  at  any  x  station  may  be  efficiently  propagated  downstream. 

Numerical  results  for  a  turbulent  wake  induced  by  a  submerged  flat 
plate  have  been  generated  and  compared  with  the  experiments  of  Swean  and 
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Keramidas.  For  pure  symmetry  boundary  conditions  for  all  variables  at  the 
free  surface,  turbulence  model  parameters  have  been  varied  so  as  to  obtain 
the  best  overall  fit  with  the  experimental  data. 

The  reference  parameter  set  seems  to  provide  a  good  description  of  the 
wake  field  for  the  cases  of  infinite  fluid.  These  are  cases  where  plate 
submergence  is  sufficient  so  that  the  free  surface  condition  does  not 
interact  with  the  flow  field  within  the  finite  downstream  domain 
considered.  For  the  finite  case,  plate  submergence  sufficiently  small  so 
that  free  surface  interaction  occurs,  the  reference  parameter  set  does  not 
adequately  describe  the  wake  field.  A  better  description  for  this  finite 
case  is  provided  by  the  parameter  set 


c 

-  0.09. 

u 

0 

O 

• 

CM 

B 

e 

% 

=  0.75, 

Cle 

-  1.43, 

C2e 

=  1.92. 

This  set  differs  from  the  reference  set  only  with  respect  to  the  effective 

Prandtl  numbers  a  and  a.  which  model  the  diffusion  terms  in  the  k 
e  k 

and  e  equations.  Diffusivity  of  e  is  retarded  while  diffusivity  of  k  is 
enhanced  relative  to  the  reference  set.  These  changes  which  are  necessary 
only  in  the  finite  case  would  appear  to  be  associated  with  free  surface 
proximity  effects . 

The  special  boundary  condition  for  e  at  the  free  surface  (3.3)  which 
has  been  found  to  be  desirable  for  applications  to  open  channel  flow  was 
introduced.  This  enhanced  e  boundary  condition  is  necessary  to  drive  the 
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effective  eddy  viscosity  to  zero  near  the  free  surface  as  observations 
predict  for  open  channel  flow.  To  obtain  agreement  with  the  experiments  on 
wake  flows  this  boundary  condition  has  to  reduce  e  at  the  free  surface  in 
opposition  to  what  is  required  for  open  channel  flows.  This  reduced  e  at 
the  free  surface  tends  to  increase  the  eddy  viscosity  at  the  free  surface. 
Applications  of  the  e  boundary  condition  which  reduce  eddy  viscosity  at  the 
free  surface  deteriorate  the  agreement  between  numerical  predictions  and 
experimental  velocity  profiles. 

As  a  consequence  of  this  behavior  for  the  imposed  e  boundary  condi¬ 
tions,  surface  proximity  effects  within  the  algebraic  stress  model  were  not 
analyzed  any  further.  Their  effect  can  be  shown  to  be  similar  to 
increasing  e  at  the  free  surface. 

From  the  results  described  above,  we  suggest  that  a  best  compromise  for 
a  turbulence  model  for  wake  decay  the  vicinity  of  a  free  surface  to  be  a 
standard  k  -  e  model  with  a  modified  parameter  set  as  given  above  and 
symmetric  boundary  conditions  for  all  variables  at  the  free  surface.  This 
is  in  agreement  with  the  proposition  of  Skop  (1).  The  model  does  not 
predict  the  experimentally  observed  migration  of  the  wake  features  toward 
the  free  surface  within  the  domain  of  the  experiment  for  the  finite  fluid 
case.  More  definitive  experiments  are  needed  to  resolve  this  behavior. 

This  is  particularly  important  since  the  predicted  velocity  profiles  with 
reference  parameters  provide  excellent  agreement  with  the  observed  profiles 
if  the  predictions  are  shifted  vertically  in  accordance  with  the  measured 
migration. 
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Figure  1  Geometry  of  the  submerged  flat  plate. 


From  Turbulence  Structure  and  Transport  Mechanism  at  The  Free  Surface  in  an 
Open  Channel  Flow,  figure  3  on  p.  515  from  the  International  Journal  of  Heat 
and  Mass  Transfer  Vol.  25,  No.  4  by  Satoru  Komori,  Hiromasa  Ueda  and 
Fumimaru  Ogino  and  Tokuro  Mizushina;  published  by  Pergamon  Press  Copyright 
©  1982;  used  by  permission. 

Figure  2  The  budget  of  the  turbulence  kinetic  energy  from  the  experiments 

of  Komori  et  al  (5). 
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From  Turbulence  Models  and  Their  Application  in  Hydraulics,  figure  16  p.  56 
by  Wolfgang  Rodi,  published  by  International  Association  for  Hydraulic  Research 
Copyright  ©  1980;  used  by  permission. 

Figure  3  Velocity  profiles  in  a  plane  wake  behind  a  flat  plate* 
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infinite  fluid  case. 
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Keraraidas  (25)  for  the  infinite  fluid  case 
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Figure  6  Turbulent  shear  stress  profiles  of  Swean  and  Keramidas  (25)  for 

the  infinite  fluid  ease. 
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Figure  7  Measured  velocity  profiles  of  Swean  and  Keramidas  (25)  for  the 

finite  fluid  case. 
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Figure  8  Measured  turbulence  kinetic  energy  profiles  of  Swean  and 

Keraraidas  (25)  for  the  finite  fluid  case. 


(mo)  A 


46 


Figure  9  Measured  shear  stress  profiles  of  Swean  and  Keramidas  (25)  for 

the  finite  fluid  case. 
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Figure  10  Predicted  velocity  profiles  for  infinite  fluid  with  reference 
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Figure  11  Predicted  velocity  profiles  for  displaced  infinite  fluid  with 

reference  parameters. 
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Figure  12  Predicted  velocity  profiles  for  Infinite  fluid  with 

o.  =  .75  o  =  2.0. 
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Figure  13  Predicted  velocity  profiles  for  displaced  infinite  fluid  with 
o.  =  .75  o  =  2.0. 
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Figure  14  Turbulence  kinetic  energy  profiles  for  infinite  fluid  with 
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Figure  15  Turbulent  shear  stress  profiles  for  infinite  fluid  with 
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Figure  16  Turbulence  kinetic  energy  profiles  for  infinite  fluid  with 
o.  =  .75  a  =  2.0. 
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Figure  17  Turbulent  shear  stress  profiles  for  infinite  fluid  with 
o,  =  .75  o  =2.0. 
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Figure  18  Velocity  profiles  for  finite  fluid  with  reference  parameters  and 

symmetric  free  surface  boundary  conditions. 
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Figure  19  Velocity  profiles  for  finite  fluid  witho^  =  .75  o£«=  2.0  and 
symmetric  free  surface  boundary  conditions. 
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Figure  20  Velocity  profiles  for  finite  fluid  with  o,  =  1.0  a  «  2.0  and 

k  e 

symmetric  free  surface  boundary  conditions. 
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Figure  21  Velocity  profiles  for  finite  fluid  witho^  =  1.0  o£  =  1.0  and 
symmetric  free  surface  boundary  conditions. 
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Figure  22  Velocity  profiles  for  finite  fluid  wltho^  =  .75  a  =  1.287  and 
symmetric  free  surface  boundary  conditions. 
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Figure  23  Velocity  profiles  for  finite  fluid  witho,  =  1.287  a  -  1.287 

k  e 

and  symmetric  free  surface  boundary  conditions. 
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Figure  24  Velocity  profiles  for  finite  fluid  witho^  =  1.0  o^  =  1.5  and 
symmetric  free  surface  boundary  conditions. 
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Figure  25  Velocity  profile  for  finite  fluid  with  a ^  =  1.0  0£  =  1.75  and 
symmetric  free  surface  boundary  conditions. 
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Figure  26.  Velocity  profiles  for  finite  fluid  with  C  =  1.43  C.  =  1.75 

16  it 

and  symmetric  free  surface  boundary  conditions. 
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Figure  27  Velocity  profiles  for  finite  fluid  with  C  =  1.43  C„  =  2.25 
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Figure  28  Velocity  profiles  for  finite  fluid  with  C,  =  1.0  =  1.92 

le  2e: 

and  symmetric  free  surface  boundary  conditions. 
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Figure  29  Velocity  profiles  for  finite  fluid  with  C.  =  1.75  C_  =  1.92 

l£ 

and  symmetric  free  surface  boundary  conditions. 
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Figure  30  Turbulence  kinetic  energy  profiles  for  finite  fluid  with 
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Figure  31  Turbulent  shear  stress  profiles  for  finite  fluid  with  reference 
parameters  and  symmetric  free  surface  boundary  conditions. 
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Figure  32  Turbulence  kinetic  energy  profiles  for  finite  fluid  with 
ok  *  .75  °L  53  2.0  and  symmetric  free  surface  boundary 
conditions. 
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Figure  33  Turbulent  shear  stress  profiles  for  finite  fluid  with 
o k  *  .75  =  2.0  and  symmetric  free  surface  boundary 

conditions. 
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Figure  38  Finite  fluid  velocity  profiles  for 
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Figure  39  Finite  fluid  velocity  profiles  for 
A  =  0.1  o  -  .75  a  =  2.0 
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Figure  AO  Finite  fluid  velocity  profiles  for 
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Figure  41  Finite  fluid  velocity  profiles  for 
A  =  1.5  o.  =  .75  o  =  2.0 
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Figure  42  Finite  fluid  velocity  rpfiles  for 
A  =»  5.0  a  =  .75  a  =  2.0 
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APPENDIX  A.  Modeling  of  the  Reynolds  Stress  Equation 

Before  we  proceed  with  the  generation  of  an  algebraic  stress  model  from 
the  Reynolds  stress  equation  certain  terms  in  that  equation  need  to  be 
modeled.  Rewriting  equation  (2.6)  in  symbolic  form  with  the  viscous 
diffusion  term  neglected  as  indicated  earlier  we  have 


where  the  stress  production  term  P^j  is  given  by 


3U,  3 


3  U 


i 


u±\  jzr  ~  ui\  37 

k  J 


k 
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large  Reynolds  numbers  anticipated  the  turbulence  is  expected  to  be 
isotropic  (11).  Therefore  is  usually  modeled  as 


ij 


3  £  6ij* 


(A.  2) 


The  pressure  strain  term  ir^  has  been  shown  by  Rotta  (17)  and 
Lumley  (12)  to  contribute  to  two  processes  one  involving  only  fluctuating 
velocities  and  the  other  involving  the  interaction  of  fluctuating  velocities 
with  the  mean  strain.  These  are  denoted  by  it  .  and  it  _  respectively. 

ij  y  1  1  j  ,  2. 

Following  Rotta  (17)  tt  is  usually  taken  to  be  proportional  to  the  local 

y 

anisotropy  of  the  turbulence.  That  is 


"ij.i  ■ '  ci  ir  (uiuj  -!kSi3) 


(A.  3) 


For  it  we  take  the  model  of  Launder,  Reece,  and  Rodi  (18) 
i  j  y  2 


9  U.  9  U , 
Y  + 

j  i 


6  (V3PV 


(A.4) 


where  D^s  -  (u^  ^  +  Ujuk  ^ 


9  U,  _  9  U,  . 

+  u,u,_  *—— )  and  P  =  y  • 
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Launder,  et  al  (18)  further  found  that  in  general  the  first  term  above 
is  dominant  and  therefore  recommended 


ij  9  2 


-  a 


( v  i 


V 


(A. 5) 


which  relates  tt^  0  to  the  anisotropy  of  the  turbulence  production  j  . 
Further,  they  introduced  wall  proximity  effects  into  the  pressure-strain 
model  by  taking  the  model  coefficients  to  be  functions  of  distance  from  the 
walls.  Daly  and  Harlow  (19)  and  Shir  (8)  introduced  an  additional  proximity 
term,  ir^  g  to  the  pressure  strain,  tt  ^  .  Following  Shir  (8)  we  take 


TT 


ij,s 


Vi  Vj  ■  f  Vj  Vi)  1 


(A. 6) 


where  f  is  a  surface  damping  function  which  decays  with  increasing  distance 

from  the  surface  and  n^  is  a  unit  normal  to  the  surface.  The  effect 

of  ii,,  is  seen  to  reduce  the  normal  fluctuations  relative  to  those 
ij  ,s 

transverse  while  conserving  turbulent  kinetic  energy.  Conservation  of  k  is 

guaranteed  since  the  trace  of  tt  vanishes.  Therefore  in  the  vicinity  of 

ij  >s 

the  surface  tt  _  retards  the  normal  turbulent  fluctuations  and  distributes 
ij  ,s 

them  adiabatically  among  the  transverse  fluctuations.  This  represents  the 
free  surface  damping  effect  eluded  to  earlier  which  has  been  observed 


experimentally  in  the  vicinity  of  a  free  surface.  Therefore,  tt  has 

11  ,  s 
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recently  been  employed  to  model  free  surface  proximity  by  Naot  and  Rodi  (6) 
and  Celik,  Hossain,  and  Rodi  (9)  in  applications  to  open  channel  flows. 

We  will  not  describe  models  for  the  other  terms  of  the  Reynolds  stress 
equation  as  they  will  not  be  needed  in  the  remainder  of  the  work  to 
follow.  Having  incorporated  the  major  physics  associated  with  near  free 
surface  flows  in  the  modeling  of  the  pressure  strain  term  we  now  generate 
the  desired  algebraic  stress  model. 


APPENDIX  B.  The  Algebraic  Stress  Model 

If  we  examine  the  Reynolds  stress  equation  (2.6)  we  note  that 
derivatives  of  the  dependent  variables  appear  in  the  time  rate  of  change, 
convective  and  diffusion  terms.  If  these  terms  could  be  expressed  as 


algebraic  forms  of  the  dependent  variable  u^uj  an(*  lower  order  correlations, 
the  differential  equations  would  become  algebraic  and  therefore  the  local 
behavior  of  the  Reynolds  stress  could  be  approximately  determined  without 
the  computational  load  of  propagating  a  modeled  differential  equation. 
Various  procedures  have  been  utilized  to  generate  algebraic  stress  models. 

Of  these  the  most  general  one  is  due  to  Rodi  (20).  Rodi  assumes  that  the 
sum  of  the  transport  terms  (rate  of  change,  convection  and  diffusion) 


for  u^uj  is  proportional  to  that  for  the  turbulent  kinetic  energy,  k.  We 

note  that  the  sum  of  the  transport  terms  for  k  is  given  from  equation  (2.9) 

as  P-e .  Therefore  we  have  from  equation  (A.l)  with  the  inclusion  of  the 

surface  proximity  term,  tt  from  equation  (A. 6) 

I  j  y  S 


Pij  +  "ij  +  "ij.s  “  eij 


uiu1 


(P-e) 


(B.l) 
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where  the  proportionality  factor  has  been  taken  to  be  u^u^/k.  This  equation 
is  seen  to  yield  an  identity  if  we  take  its  trace  and  Rodi  (21)  suggests 


that  it  represents  a  good  approximation  when  the  variation  of  u^u^/k  is 


small  compared  to  the  variation  of  u^uj •  Equation  (B.l)  provides  the 
desired  algebraic  stress  model  after  substitution  of  the  modeled  forms 
of  ir^  ,  g  and  from  equations  (A.  3),  (A. 4)  or  (A.  5),  (A.  2)  and 
(A. 6). 

Before  we  proceed  with  these  substitutions  a  few  comments  are  in  order 
with  respect  to  the  general  algebraic  stress  equation  above.  For  flows  that 
evolve  very  slowly  it  is  often  justifiable  to  neglect  the  time  rate  of 


change  and  transport  of  u^uj •  This  then  is  equivalent  to  setting  P  ae  in 

equation  (B.l).  Sometimes  this  is  done  for  the  off-diagonal  stress 

components  only  with  the  normal  components  evaluated  with  the  full  algebraic 

equation  P  ^  £  (9)  because  of  their  close  relationship  to  k. 

When  convection  and  diffusion  are  neglected  completely,  corresponding 

to  P=e  in  equation  (B.l),  and  after  substitution  of  the  modeled  terms 

for  7T  ,  7T  and  e.,  from  equations  (A. 2),  (A. 4)  and  (A. 6),  the  following 

ij  ij 

equation  results 


a  -  a)  PtJ  -6  Dy-T  k(5^-  +  8J^  *!'(■  +*>«„ 


V! 


(B.  2) 


+  C!  r-  (u2  6  -  u  u.6  -  1  u  u.  6  )  f  =  0 
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where  as  defined  previously, 


The  last  term  of  equation  (B.2)  is  recognized  as  Shir's  (8)  model  applied  to 
the  free-surface  proximity.  The  variable  in  that  term  represents  the 
fluctuating  velocity  component  normal  to  the  free  surface. 

The  algebraic  stress  relation  given  above  has  been  developed  for  a 
general  turbulent  flow  field  in  the  vicinity  of  a  free  surface.  Further 
simplifications  can  be  made  if  we  introduce  some  of  the  specialized  features 
of  the  wake  application  to  be  pursued  here.  First,  we  anticipate  the  flow 
field  to  be  predominantly  in  the  x  direction,  so  that  the  velocity 
components  in  the  y  and  z  direction  are  expected  to  be  much  smaller  and  in 
fact  of  secondary  nature  (unidirectional  flow).  In  addition,  with  respect 
to  a  ship-fixed  coordinate  system,  the  basic  flow  is  expected  to  be 
steady.  Then  the  asymptotic  analysis  of  Skop  (1)  shows  for  large  Reynolds 
number  that  diffusion  in  the  x  direction  will  be  negligible  compared  to  that 
in  the  y  and  z  directions.  That  is  to  say,  we  anticipate  a  parabolic  flow- 
field  description  will  be  valid.  In  addition,  this  characteristic  renders 
the  k  and  £  equations  (2.8)  and  (2.9)  parabolic  also.  This  is  certainly  to 
be  desired  so  that  a  marching  technique  can  be  applied  to  generate  the 
solution  and  no  iteration  will  be  required. 

If  we  introduce  the  assumption  of  unidirectional  flow,  U^=  U  (y,z)  e^ 
where  y  is  perpendicular  to  the  free  surface  into  the  stress  relations,  we 
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obtain  the  following  algebraic  modeling  equations  for  the  various  Reynolds 
stress  components  (16): 


1(V  2CJf)  ? 


+  6  [uv 


8U  i 

57  J 


(B.  3) 


i  ciw2  "  cif  7  ^  +  t  (a  - p  +  4  ( cr0  e 


Q  r —  ou  —  c?  ui 

“  P  [uv  r - UW  r— J 

3  y  3  zJ 


(B.4) 


Ic  ^  \  ^  uv  =*  -  [  ( 1-a  )  (  +  vw  tan  ij>) 
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—  C  _uw  =  -  [  (  1-a  )(  w^  +  vw  cot^j  -  8  1/  +Tkl  . 
k  1  L  v  J  3z 


(B.6) 


and 


(B.7) 


=  (i»)  /  (*» 

l3zJ  '  l3y 

dimensionality  of  the  basic  flow.  That  is,  for  two-dimensional  flow 


where  tan  =  (^-J  /  and  therefore  tan  characterizes  somewhat  the 
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where  ■  0  ,  tan  “  0.  The  two-dimensional  plane  wake  which  is  of 
interest  is  considered  below. 

Naot  and  Rodi  (16)  developed  an  iterative  procedure  to  solve  these 
equations  for  the  stress  components  for  fixed  values  of  the  k-e  parameters 
and  tan  and  they  have  evaluated  the  variation  of  the  stress  components 
with  the  surface  proximity  CJf  and  the  inclination  angle  ^ .  For  all 
inclination  angles  they  demonstrate  that  these  equations  provide  for  the 
damping  of  the  normal  fluctuations  and  the  enhancement  of  the  transverse 
components  as  we  approach  the  free  surface.  Additionally,  the  Reynolds 
shear  stress  component  uv  is  shown  to  be  damped  as  the  free  surface  is 
approached.  For  computational  purposes,  they  suggest  the  following 
approximations  for  the  turbulent  shear  stresses  accurate  to  ten  percent: 
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and 


-  uw  =  C 
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(B.  9) 


where  =  0.09.  These  are  in  the  form  of  eddy  viscosity  relations  with  the 
effective  eddy  viscosity  a  function  of  free  surface  proximity.  Note  that  at 
remote  locations  from  the  free  surface  f  -*■  0  and  we  recover  the  standard 
eddy  viscosity  relations 
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for  both  stress  components-  As  we  approach  the  free  surface,  f  increases 
and  the  effective  eddy  viscosity  decreases.  Different  eddy  viscosities  are 
predicted  for  the  vertical  and  transverse  transport  of  momentum  uv  and  uw 
because  of  the  surface  proximity  effect. 

The  two  equations  (B.8)  and  (B.9)  for  uv  and  uw  are  complemented  by 
equations  (B.3)  (B.4)  and  (B.7)  for  v* ,  ^  and  vw  to  form  an  algebraic 
stress  model  including  free  surface  effects.  However,  the  surface  proximity 
function  f  still  must  be  defined. 

The  surface  proximity  function  is  taken  to  be  dependent  upon  (£/y), 
where  &  is  the  local  dissipation  length  scale  of  the  turbulent  motion  at  the 
point  considered  and  y  is  the  distance  of  this  point  from  the  free 
surface.  Most  authors  have  assumed  linear  variations  of  f  with  i  /y,  e.g. 
Launder  (22),  Gibson  and  Launder  (23)  and  Celik  et  al  (9),  in  applications 
to  both  wall  proximity  and  free-surface  proximity  for  open-channel  flows. 
Naot  and  Rodi  (6)  have  taken  quadratic  variations  for  f(£/y).  A  general 
form  for  f  can  therefore  be  written  as 


where  n  =»  1,2  depending  upon  whether  a  linear  or  quadratic  dependence  is  to 
be  utilized  and  is  a  constant  which  is  to  be  specified. 

In  the  open-channel  applications,  free-surface  proximity  effects  have 
been  included  in  addition  to  simultaneous  wall  proximity  effects.  In  our 
present  application  to  the  far  wake  problem,  the  walls  which  may  have  been 
responsible  for  generating  the  wake  are  so  remote  from  the  region  of 
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interest  that  they  need  not  be  included-  Our  application  therefore 
considers  the  free  surface  proximity  only.  Once  a  choice  has  been  made  for 
n  in  equation  (B.10),  the  surface  proximity  effect  is  characterized  by  two 
constants  Cj  and  Cf. 

1 .  Two-Dimensional  Flow 

The  algebraic  stress  model  equations  given  above  have  been  developed 
for  uni -directional  three-dimensional  flow.  For  a  two-dimensional,  uni¬ 
directional  flow  such  as  the  NRL  flat  plate  experiments  (25),  the  equations 
are  simplified  as  are  the  computational  demands  to  generate  turbulent  flow 
field  predictions.  The  geometry  is  illustrated  in  Fig.  1. 

For  two-dimensional  planar  flow  we  can  neglect  all  variations  in  the  z 

g 

direction.  That  is,  r —  *  0  and  also  the  Reynolds  shear  stress  components 

d  z 

uw  and  vw  are  zero.  We  must,  however,  include  the  normal  stress  component 
w2  since  this  will  not  be  zero.  The  required  algebraic  stress  relations 
therefore  become 
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and 


§  Cjf)  uv 
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We  recognize  that  the  last  terms  in  equations  (B.ll)  and  (B.12)  are 
proportional  to  the  turbulence  production  P  for  two-dimensional  uni¬ 
directional  flow, 

p  -  W  !y.  (B.  14) 

Substituting  (B . 14 )  into  (B.ll)  and  (B.12)  and  making  the  local  equilibrium 
assumption  P  =  e  and  2k  a  u2  +  v2  +  ,the  following  solution  is  obtained: 
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and 
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(B. 17) 


For  the  case  8  =*  y  =  0  often  employed  in  these  models  as  a  consequence 
of  utilizing  the  dominant  behavior  for  ir.  -  given  in  equation  (3.5),  these 

*.1  9  & 

results  simplify  further  to  yield 
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Equation  (B.18)  is  seen  to  be  very  similar  to  the  approximation  of  Naot  and 
Rodi  (16)  given  previously  by  equation  (B.8).  It  represents  an  eddy 
viscosity  formulation  for  the  Reynolds  shear  stress  with  the  surface 
proximity  behaving  identically  to  that  of  equation  (B.8).  The  surface 
damping  of  the  normal  stress  v^  and  the  enhancement  of  is  apparent  from 
equations  (B.19)  and  (B.20).  We  can  make  a  complete  connection  of  Equation 
(B.18)  with  Equation  (B.8)  by  defining  in  (B.18) 
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2  (1  -  a) (a  +  Cj-  1) 


=  C 


This  yields 
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Note  that  in  Equation  (B.21)  in  the  limit  f  +  0  we  recover  the  standard  eddy 

viscosity  formulation.  That  is,  Eq.  (B.21)  behaves  as  a  standard  eddy 

viscosity  in  the  absence  of  surface  proximity  effects.  The  surface 

proximity  effects  are  seen  to  reduce  the  eddy  viscosity  as  the  free  surface 

is  approached.  Equation  (B.21)  is  seen  to  be  sufficient  to  close  the 

turbulence  problem  for  the  two-dimensional  case.  In  its  given  form,  it 

C1 

provides  for  two  additional  parameters  and  Cf.  Cf  is  necessary  to 

C1 

define  the  surface  proximity  function  f.  Therefore  combining  this  with  the 
k  and  e  equations  (2.8)  and  (2.9)  which  were  seen  to  contain  five  constants, 
our  complete  turbulence  model  including  free-surface  proximity  effects 
possesses  seven  constants  or  adjustable  parameters. 
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